The single-atom box: bosonic staircase and effects of parity 
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We have developed a theory of a Josephson junction formed by two tunnel-coupled Bose-Einstein 
condensates in a double-well potential in the regime of strong atom-atom interaction for an arbitrary 
total number N of bosons in the condensates. The tunnel resonances in the junction are shown to 
be periodically spaced by the interaction energy, forming a single-atom staircase sensitive to the 
parity of TV even for large N. One of the manifestations of the staircase structure is the periodic 
modulation with the bias energy of the visibility of the interference pattern in lattices of junctions. 
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Bose-Einstein condensates (BECs) have a unique abil- 
ity to give rise to individual quantum states occupied 
with a macroscopically large number of particles. The 
interference of the wavefunctions of such macroscopic 
states leads to a variety of macroscopic quantum phe- 
nomena, the most studied of which is the Josephson ef- 
fect: the flow of current across a potential barrier due to 
the interference of the condensate states on the two sides 
of the barrier. For a long time since the original discovery 
[H, this effect has been studied in solid-state "Josephson 
junctions" Qj: two weakly-coupled superconductors, in 
which the condensates are formed by Cooper pairs 
At energies smaller than the binding energy of electrons 
in Cooper pairs, the characteristics of Josephson junc- 
tions can be discussed directly in terms of the properties 
of the Cooper-pair condensates 0] ■ In this case, the two 
dynamic variables of the Josephson junction (JJ) are the 
phase difference <fi between the condensate wavefunctions 
on the two sides of the barrier and the half-difference n of 
the number of Cooper pairs in the two condensates. Since 
the phase </> and the "charge" n are conjugate quantum 
variables [5| , the dominant qualitative feature of the J J 
dynamics is the transition from the classical dynamics 
of <j) in large junctions to the classical dynamics of n in 
small junctions. The classical charge regime is character- 
ized by the correlated transfer of individual Cooper pairs 
Q, with the simplest "Cooper-pair box" system 0, [1] 
providing the most basic demonstration of this transfer 
in the "Cooper-pair staircase" : the step- wise increase of 
n with the junction bias that corresponds to successive 
transfers of Cooper pairs one by one through the junc- 
tion. The transition between the quasiclassical dynamics 
of <j) and n is affected by the energy dissipation in the 
junction 0] (with the classical-n regime requiring small 
dissipation) but can generally be understood as compe- 
tition between the tunnel amplitude and the interaction 
energy of Cooper pairs. 



In the case of solid-state JJs, the bosonic junction 
model is partly a pedagogical simplification. With ul- 
tracold atoms in optical traps of double-well structure, 
this model, however, can id. Ill| and recently has been 



realized directly. It has been studied experimentally both 
in the regime of larger-scale BECs in the wells with rel- 
atively weak transverse confinement [l2|, [l3|, when 
they exhibit the Josephson effect governed by the classi- 
cal dynamics of the phase difference <j>, and in the "Mott 
insulator" regime [HI, [3] , when the total number N of 
atoms in the double- well structure is small N = 1, 2 and 
the number difference n behaves quasiclassically. One of 
the main dynamic features of the regime of classical n 
is the existence of resonances, as a function of the en- 
ergy bias e, cor resp onding to the tunneling of individual 



atoms 



The goal of this work is to study 
resonances for an arbitrary number 

in both 



these "single- atom" 

of atoms N in the BEC Josephson junction, 
the interaction-dominated Mott-insulator regime and the 
case of large-scale BECs. We show that the main features 
of the single-atom effects in BEC junctions are quite sim- 
ilar to those in solid-state JJs [H, 0, HI , and in particular, 
include the periodic spacing of resonances in the "single- 
atom staircase" . There are, however, several new effects 
produced by the fixed total number of bosons N in the 
BEC junction, as opposed to solid-state JJs where the 
total number of particles is not well-defined. One is the 
modulation of the strength of the resonances with the 
number difference n, which reflects the collective nature 
of the tunnel coupling of the BECs even in the regime 
of single-atom tunneling. Another is the dependence of 
the resonances on the parity (even or odd) of the total 
number of bosons N. 



We consider N bosons in a double-well potential ap- 
proximated by a two-site Bose-Hubbard model [13, 21 
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22, 22 
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H = — ^2 ni(n l -l)-A(ala2+a2ai)~e(n 1 ~-n 2 ). (1) 

i=l,2 

Here rii — a\di and a\ {a/) creates (annihilates) a bo- 
son in the lowest-energy state localized in the ith well, 
A is the tunnel coupling of these two states, 2e is their 
difference in on-site energy, and U is the on-site interac- 
tion. The Hamiltonian fl} conserves the total number of 
bosons N — m + n 2 , and after separation of N, it takes 
the following form (up to a constant energy shift): 



Hq = Un 2 — 2en — A(a\a,2 + a\a\) . 



(2) 



Here n = (m — U2)/2 is the occupation number differ- 
ence between the two wells. While N does not enter 
the Hamiltonian |J2) explicitly, it influences the junction 
dynamics by limiting the range of the possible values of 
the number difference n € {-N/2, -N/2 + l,...,N/2}, 
i.e., to integer or half- integer n for even or odd N, re- 
spectively. In general, the kinematics of the model de- 
fined by Eq. {1} or Q gives Schwinger's representation 
of the angular momentum J of magnitude N/2 (see, e.g., 
[2~ij]). In this identification, n is the z-component J z of 
the momentum, and the tunneling term in the Hamilto- 
nian is —2AJ X . Since the properties of the integer and 
half-integer momenta are quite different, this parity effect 
makes the energy characteristics of the bosonic junction 
([2|) dependent on the parity of N, as we show in more 
detail below. 

To see this quantitatively, we start with the regime of 
strong interaction, U 3> NA. In this case, the tunnel 
coupling A = A[(N/2 + n+ l){N/2 - n)] 1 / 2 of the states 
n and n + 1 that differ by the transfer of one atom in the 
junction, is small, A -C U. This means that these states 
are coupled effectively only in the vicinity of certain val- 
ues of the energy bias e s» e„, 



(n+l/2)U, 



(3) 



when the interaction energy difference between them is 
suppressed and atoms can tunnel across the junction. 
Since these resonances are uniformly spaced by U, the 
average number difference (n) in the ground state of the 
junction shows a regular staircase-like structure (Fig. [T]) 
similar to the corresponding structure in the Cooper-pair 
box @, The curves in Fig. [T] are obtained by straight- 
forward numerical diagonalization of the Hamiltonian ([2]) 
and show also the smearing of the staircase by tunneling 
when it is increased to the level A ~ U/N. At U > NA, 
the system dynamics reduces in the vicinity of each step 
to the two states n and n + 1, and one can find directly 
the shape of the steps in the staircase: 



(n) 



[(e-e„) 2 + A 2 ]V2 



(4) 
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FIG. 1: Average number difference (ri) in the ground state 
of the double-well BEC Josephson junction for N — 3 as a 
function of the energy difference e between the two wells, for 
several values of the ratio of interaction energy versus tun- 
nel amplitude U/A. Curves with U/(NA) > 1 show the 
single-atom staircase, which is smeared with increasing tun- 
neling strength. Left inset: Compressibility n = d(n)/de 
showing the peak structure associated with the single-atom 
resonances. For odd N, there is a resonance in the unbi- 
ased junction e ~ 0. Right inset: k for = 10 bosons in 
the two wells. With increasing tunneling strength, k, changes 
from sharp single-atom resonances to a quasiclassical behav- 
ior with constant large-scale compressibility ~ 1/(7 and weak 
sinusoidal single-atom oscillations. 



The shape of each of the single-atom transitions is seen 
more clearly in the resonant peaks in the compressibility 
K = d(n)/de of the ground state. In particular, Eq. ((4]) 
shows that the maximum of each peak reflects directly 
the effective tunnel coupling of the two states: 



1 

2A 



= —[(N/2 + n+l)(N/2-n)]- 1 / 2 . (5) 



This equation agrees with the curve for strongest inter- 
action in the right inset to Fig. [TJ and shows that the 
enhancement of the tunneling amplitude from the single- 
particle value A to the multi-particle value A ~ NA by 
the BEC coherence manifests itself even in the regime of 
single-atom tunneling. The modulation ([S]) is a quantum 
consequence of the dependence of the critical current of 
the BEC junction on the population imbalance between 
the two wells that also leads to "self-trapping" of the 
condensate in the case of classical Josephson dynamics 

IlGIEl. 

The main qualitative feature of the staircase that re- 
flects the parity effect is that the steps ((4]) occur at half- 
integer values of energy bias e (in units of interaction 
energy U) for even N and integer e/U for odd N, see 
Eq. ([3|) . In particular, the tunneling is always suppressed 
for even N in the symmetric unbiased junction, when 
e ~ 0, but there is a tunneling resonance at e = for odd 
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N (Fig. [J). This resonance at e = has an obvious inter- 
pretation in the single-atom regime with strong interac- 
tion energy, U 3> AA: the interaction energy of N atoms 
is degenerate with respect to the last atom out of an odd 
total number N being placed in the left or right well. 
This simple interpretation is not valid in junctions with 
large N, when U <§C NA, but the parity effect exists even 
in this case, provided the interaction energy still plays a 
role in the junction dynamics, U ^S> A/N. (If U <C A/AT, 
the interaction is not essential and the junction repre- 
sents a collection of N independent two-state systems of 
individual atoms.) In this regime, the number-difference 
fluctuations Sn around the average value (n) = e/U have 
the magnitude (estimated from the Hamiltonian ([6]) be- 
low) 1 < Sn ~ (A/A/f7) 1/4 < VW. If (n) is not too 
close to the boundaries ±N/2 of the allowed interval for 
n, this means that effectively there are no restrictions 
on the values of n, and variations of the tunneling am- 
plitude A are negligible. As a result, one can define a 
phase 4> € [0, 2n] conjugate to n and the Hamiltonian ([2|) 
takes a form that essentially coincides with that of the 
"solid-state" Josephson junction: 



H = — 



U 



u 



o 2 



2A(e) cost 



(6) 



where A(e) = A[(N/2) 2 - (e/U) 2 ] 1 / 2 . The boundary con- 
ditions on the wavefunctions ip {(/>), however, depend on 
the bias e and the parity of A": 



(7) 



Using the known result [6( for the Hamiltonian ([6]) we see 
that the ground state energy of the bosonic junction is 



E = -(e z /U) - (-l) iv <5 cos(27re/C/), 



(8) 



where S ~ (l/ 7 r)(At/) 1 / 2 exp[-8(A/J7) 1 / 2 ]. The first 
term in (JSj) describes the large-scale, mean-field com- 
pressibility k = —(\/2)d 2 Eo/dt 2 = 1/U for e within the 
range [-NU/2,NU/2] (right inset to Fig. [J). Qualita- 
tively, this compressibility is dominated by the interac- 
tion, since the total interaction energy scales as N 2 as 
opposed to the tunneling energy which is proportional to 
N. The second term in ((HJ describes the sinusoidal single- 
atom oscillations of small amplitude which increase at the 
ends of the range of e because of the decrease of A (right 
inset to Fig.[l}. These oscillations represent what is left 
of the single-atom resonances when they are washed out 
by strong tunneling, i.e., describe the discreteness of the 
BEC flow in this quasiclassical regime. The position of 
the oscillations as a function of the bias e depends on the 
parity of N, with a minimum in n at e ~ for even AT 
and a maximum for odd N. 

The sensitivity of the bosonic JJ to individual atoms 
makes it interesting to consider the situation with one ex- 
tra "foreign" atom in the junction. Motivated by exper- 
iments on Bose-Fermi mixtures [26[ , we assume that this 



extra particle is a "fermion" , although its statistics does 
not play a role in our discussion. Besides the bosonic 
part Ho ((TJ, the junction Hamiltonian includes then a 
fermionic part Hf of the same structure with the tunnel 
amplitude A/, the number operator n/j, on-site interac- 
tion Uf, and a boson- fermion interaction term: 



H = H + Hf + Ubf(nirifi + n 2 nf 2 ) 



(9) 



The boson- fermion interaction energy Ubf can be positive 
or negative [26j |. The difference e in on-site energies is 
assumed to be the same for bosons and the fermion. 

One can follow separately the dynamics of the two 
types of particles, and we limit our brief discussion to 
the behavior of the bosons. The average boson number 
difference (n) is shown in Fig. [2] for AT = 3. In the case of 
vanishing boson-fermion interaction, Ubf = 0, we obtain 
the familiar three steps of height 1. For strong attrac- 
tion, Ubf <C —U, or repulsion, Ubf ^> U, minimization of 
the interaction energy requires that the bosons are mov- 
ing all together, and (n) shows only one step of height 
3 at e = 0. In our example of 3 bosons, this case is 
realized already for \Ubf\ > 2U . In the intermediate re- 
gion \Ubf\ < 2U, we find a more complex behavior: for 
Ubf < and Ubf > the bosonic steps are shifted, re- 
spectively, to and away from the region e = 0. Moreover, 
for Ubf > (i.e., for antiferromagnetic coupling between 
the angular momenta representing bosonic and fermionic 
parts of the junction) , the system is frustrated and the 
steps have regions with negative slope, see the curves 
for Ubf jU = 0.5 and 1. In the "single-atom" regime of 
small tunneling amplitudes, and for odd N, the region 
of the negative slope around e ~ corresponds to the 
range e € [—Ubf/2, C/&//2], where (n) is described by the 
following general relation: 
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2 [e 2 (A 2 
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A 2 ) 2 + 



\ 2 A 2 U 2 fY/ 2 



(10) 



Note that the slope depends of the relation between the 
effective boson and fermion tunneling amplitudes. 

One way of studying the effects discussed in this work 
would be to measure the interference pattern of the ex- 
panding atomic cloud after switching off an optical lattice 
of appropriately spaced uncoupled double- well junctions. 
Such arrays of double wells have already been realized 
experimentally [271 ]. The interference of particles within 
each junction will produce a cosine modulation of the av- 
erage atom intensity /. The visibility of this modulation 
pattern is given [281 ] by the average kinetic (tunneling) 
energy per boson in units of A: 



{In 



l )/(In 



t ) = (a\a 2 + a\ai)/N . 



In the single-atom regime of the bosonic junction 
the visibility is modulated by the bias e and reaches its 
maximum at the resonances Q 

-1/2 



v= {\/N)[{e-e n Y + \ l 



(11) 
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FIG. 2: (color online) Steps in (n) as a function of e for N = 3 
bosons, iV/ = 1 fermions, A/Z7 = 0.01, and A f /U = 0.004. 



where the states of the two wells deviate most strongly 
from pure number states. For the Bose-Fermi junction 
([9]), v shows a rich dependence on e and Ubf that is also 
related to the staircase structure, see Fig. [3J 
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FIG. 3: (color online) Visibility v of the bosonic interference 
pattern in the Bose-Fermi junction. All the parameters are 
chosen as in Fig. [2] 

In conclusion, we have shown that a Josephson junc- 
tion formed by two tunnel-coupled BECs shows single- 
atom effects similar to those in solid-state Josephson 
junctions. In addition to periodically spaced resonances 
as a function of the difference in on-site energy leading 
to a single-atom staircase, there are modulations of the 
strength of the resonances and parity effects that can be 
traced back to the fixed total number of bosons N in 
the BEC junction. A different, e.g. fermionic, additional 
particle in the junction leads to non-trivial modifications 
of the staircase, that can be experimentally observed in 
the visibility of the interference pattern. 
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